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1. Introduction 

The theoretical concepts of Brownian motors and ratchet transport pQ have been 
experimentally realized in a variety of systems. Examples are: cold atoms in optical 
lattices f2\ . colloidal particles in holographic optical trapping patterns ratchet cellular 
automata superconducting films with periodic arrays of asymmetric pinning sites 
El, to mention only a few. 

When we study the motion of Brownian motors, the natural transport measure 
is a conveniently defined average asymptotic velocity (v) of the Brownian motors. It 
describes how much time the typical particle needs to overcome a given distance in 
the asymptotic (long-time) regime. This velocity, however, is not the only relevant 
transport criterion. Other attributes can also be important. In order to establish these, 
we consider the two following aspects: the quality of the transport and the energetic 
efficiency of such a system. 




Figure 1. (Color online) Two sets of illustrative trajectories of an inertial, rocked 
Brownian motor (see in text). Both sets of trajectories A and B possess the same 
average asymptotic velocity, but exhibit a distinct different diffusion behavior. 

In Fig. ^ one can identify two different groups A and B of random trajectories of 
the Brownian particle; both possess the same average drift velocity {v). However, it is 
obvious upon inspection that the dynamical properties of these two groups of trajectories 
are different. The particles from the group A travel more or less coherently together 
while the particles from the group B spread out as time goes by. If we fix the distance 
x = Xi then most particles from the group A reach this distance at about the same time 
t = ti, while for t = t± most of the B trajectories either stay behind or have already 
proceeded to more distant positions. It is thus evident that the noise-assisted, directed 
transport for the particles in the group A is more effective than in the group B. 
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There is still another efficiency aspect related to Brownian motor transport. This 
refers to the external energy input into the system which may be essential in practical 
applications. We like to know how much of this input energy is converted into useful 
work, namely into directed cargo transport, and how much of it gets wasted. Since 
motors move in a dissipative environment, we need to know how much of the input 
energy is being spent for moving a certain distance against the acting friction force. 
Fig. El depicts trajectories representing different motor scenarios. The motor C moves 
forward unidirectionally, while the motor D moves in a more complicated manner: its 
motion alternates small oscillations and fast episodes, mostly in the forward direction, 
but sometimes also in the backward direction. Again the mean velocity in both cases 
is the same, however, the particle C uses energy pumped from the environment to 
proceed constantly forward while the particle D wastes part of its energy to perform 
oscillations and back-turns. By simply inspecting these schematic pictures one can guess 
immediately when directed transport is more effective. 



t 






X 


C 








jT ► 


X 


D 


>■ 



t 

Figure 2. (Color online) Typical trajectories of an inertial, rocking Brownian motor; 
both sets assume the same average velocity but differing velocity fluctuations. 

We note that in Fig. ^ the cases A and B can be characterized by the effective 
diffusion coefficient -D e //, i.e., by the spreading of fluctuations in the position space 
while the cases C and D in Fig. 2 can be characterized by the variance of velocity 
a v = ( v2 ) ~ ( v ) 2 - The three quantities (v), D e ff and a% can be combined to define two 
important characteristics of transport, namely the efficiency of noise rectification and 
the so-called Peclet number 

Our work is organized as follows. In the following section, we detail the model 
of an inertial rocked Brownian motor. In section 3, we present a general discussion of 
the efficiency measures of Brownian motors. In section 4, a description of the ratchet 
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based on point processes is introduced. In sections 5 and 6, our numerical findings are 
analyzed in the context of the optimization conditions for transport of inertial Brownian 
motors. A summary is provided in section 7. 

2. Inertial rocked Brownian motors 

The archetype of the inertial Brownian motor is represented by a classical particle of 
mass m moving in a spatially periodic and asymmetric potential V(x) = V(x + L) with 
period L and barrier height AV [HI E] • The particle is driven by an external, unbiased, 
time-periodic force of amplitude A and angular frequency Q (or period To = 2ir/Q). 
The system is additionally subjected to thermal noise The dynamics of the system 



V(x,t) 




X 

Figure 3. (Color online) Schematic picture of a rocking ratchet with the potential 
V(x,t) = V(x) - xacos(ujt), cf. Eqs (2) and (4). 

is modeled by the Langevin equation [TU] 

mx + 7s = -V'(x) +Acos(Qt) + <j2~/k B T f (t), (1) 

where a dot denotes differentiation with respect to time and a prime denotes a 
differentiation with respect to the Brownian motor coordinate x. The parameter 
7 denotes the Stokes friction coefficient, ks the Boltzmann constant and T is the 
temperature. The thermal fluctuations due to the coupling of the particle with 
the environment are modeled by a zero-mean, Gaussian white noise £(t) with auto- 
correlation function (s)) = 5(t — s) satisfying Einstein's fluctuation-dissipation 
relation. 

Upon introducing characteristic length scale and time scale, Eq. (0) can be 
rewritten in dimensionless form, namely 

l + j£ = -V\x) + a cos(ai) + y^A, i(t), (2) 

with [TT] 

x - t a mL 2 
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The characteristic time r is the time a particle of mass m needs to move a distance 
L/2 under the influence of the constant force AV/L when starting with velocity zero. 
The remaining re-scaled parameters are: 

• the friction coefficient 7 = (7/771)70 = tq/tl is the ratio of the two characteristic 
times, To and the relaxation time of the velocity degree of freedom, i.e., 77 = m/7, 

• the potential V(x) = V(x)/AV = V(x + 1) has unit period and unit barrier height 
AV = 1, 

• the amplitude a = AL/AV and the frequency u = Qt (or the period T = 2tt/uj), 

• the zero-mean white noise £(£) has auto- correlation function (£(t)£(s)) = S(t — s) 
with re-scaled noise intensity D = ksT/AV. 

From now on, we will use only the dimensionless variables and shall omit the "hat" for 
all quantities in Eq. (J2J). 

For the asymmetric ratchet potential V(x) we consider a linear superposition of 
three spatial harmonics [11, 

V(x) = Vo[sin(27rx) + c\ sin(47rx) + C2 sin (67m:)], (4) 

where Vq normalizes the barrier height to unity and the parameters C\ and C2 determine 
the ratchet profile. Below, we analyze the case when c\ = 0.245 and C2 = 0.04. Then 
Vq = 0.461. This potential is shown as a bold (red) line in Fig. El 



3. Quantifiers characterizing optimal transport of Brownian motors 

As already elucidated above, there are several quantities that characterize the 
effectiveness of directed transport. The effective diffusion coefficient, describing the 
fluctuations around the average position of the particles, is defined as 

g „,= ft- -<*<«»', (5) 

" t^co 2t w 

where the brackets (...) denote an average over the initial conditions of position and 
velocity and over all realizations of the thermal noise. The coefficient D e ff can also 
be introduced via a generalized Green-Kubo relation which we detail in the Appendix. 
Intuitively, if the stationary velocity is large and the spread of trajectories is small, the 
diffusion coefficient is small and the transport is more effective. To quantify this, we 
can introduce the dimensionless Peclet number Pe ^] by use of a double-averaging 
procedure, i.e., 

V eff 

where the "double-average" ((v)) denotes the average of the asymptotic velocity over 
one cycle of the external drive, i.e., 

((«)> = lim \ f\v(t')) alt* = r'\v(t')) as dt', (7) 



t^oo t Jo 2tx Jo 

where the average (. . .) as in the second integral refers to the asymptotic periodic state. 
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Originally, the Peclet number Pe arises in problems of heat transfer in fluids where it 
stands for the ratio of heat advection to diffusion. When the Peclet number is small, the 
random motion dominates; when it is large, the ordered and regular motion dominates. 
The value of the Peclet number depends on some characteristic length scale of the 
system. Dealing with ratchets the most adequate choice for such length scale is the 
period of the periodic potential, which in re-scaled units is equal to 1. 

The second aspect of the motor trajectories we want to control has to do with the 
fluctuations of the velocity v(t). In the long-time regime, it is characterized by the 
variance o~\ = ((v 2 )) — ((f)) 2 . The Brownian motor moves with an actual velocity v(t), 
which is typically contained within the interval 

v(t)E(((v))-a v ,((v)) + a v ). (8) 

Now, if cr„ > ((f)), the Brownian motor may possibly move for some time in the direction 
opposite to its average velocity {{v)) and the directed transport becomes less efficient. If 
we want to optimize the effectiveness of the motor motion we must introduce a measure 
for the efficiency rj that accounts for the velocity fluctuations, too, namely ^3] 

V \((vy + o*-D \ \((v 2 ))-D \- {) 
This definition follows from an energy balance of the underlying equation of motion (J2J) 
(see the Appendix of Ref. [TT]). If the variance of velocity a v is reduced, the energetic 
efficiency © increases and the transport of the Brownian motor becomes more efficient. 

4. A corresponding point process related to the rocked Brownian motor 
dynamics 

The running trajectories can be characterized in a coarse grained way by only counting 
the events when a trajectory traverses from one potential well into a neighboring one, 
and by disregarding the details of the intra-well motion, see figure 0] In this way, 
a point process can be introduced that can be investigated in a standard way [T4*j . 
Most, though not all, of the quantities describing the original continuous process can be 
retrieved from the so defined point process. For this purpose we introduce two random, 
natural numbers N£, where a = {R,L} stands for right (R) and left (L). The number 
Nj* is given by the number of barrier crossings towards the right within the fc-th period 
of the driving, i.e. in the time between (k — 1)T and kT . The respective number of 
barrier crossings to the left is denoted by N%. The difference 

N k = JV* - N% (10) 

indicates that during a temporal period T the particle has covered the distance 

Xk = AfcL = Nk- Hence the average, asymptotic velocity is given by 

i r t 1 K r kT 

M = lim - £ WO) <* = Jim — g jf^MO) * 

- £ t - ±- { I n„ = M. (id 
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Figure 4. (Color online) The point process related to rocked Brownian motor 
dynamics. Two thresholds LT and RT are the corresponding maxima, to the left and 
to the right of the particle position. If the particle jumps over the right neighboring 
barrier, the number increases to Nj^ + 1 and the previous right sided threshold 
assumes the role of the left sided threshold. 



Analogously, the effective diffusion coefficient is determined by the relation 

m= OT-w» (12) 

A related quantity is the Fano factor F ^S], defined here as the fluctuation to the first 
moment ratio 

As such, the Fano factor provides a quantitative measure of the relative number 
fluctuations or the relative randomness of the process; in the case of a Poisson process 
F = 1. 

On the other hand, from (JSj). (jllj) and (fT2|) it follows that the Peclet number can 
be expressed as 

Pe=^- (14) 
This quantifier is thus related to the Fano factor via the relation Pe = 2/F. 

5. Numerical analysis 

The noiseless, deterministic inertial rocked ratchet shows a rather complex behavior 
and, in distinct contrast to overdamped rocked Brownian motors ^H], often exhibits a 
chaotic dynamics [HI EI- By adding noise, one typically activates a diffusive dynamics; 
thus allowing for stochastic escape events among possibly coexisting attractors. As 
analytical mathods to handle these situations effectively do not exist, we carried out 
extensive numerical simulations. We have numerically integrated Eq. (J2J) by the Euler 
method with time step h = 5 x 10~ 4 T. The initial conditions for the coordinate x(t) 
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were chosen according to a uniform distribution within one cell of the ratchet potential. 
The starting velocities of the particles were also distributed uniformly in the interval 
[-0.2,0.2]. 

The first 10 3 periods T of the external force were skipped in order to avoid 
transient effects. We employed two tactics of extracting the above characteristics from 
the generated trajectories. For the estimation of the energetic efficiency (or velocity 
fluctuations) the usual averages over the time (10 5 T) and 333 different realizations 
were taken. In the the case of the Peclet number (or effective diffusion) we used the 
point process approach (see previous section for details); therefore, only one time average 
of long-time runs (10 6 T) was required. 

Typically, there are two possible dynamical states of the ratchet system: a locked 
state, in which the particle oscillates mostly within one potential well (cf. the case with 
amplitude a = 4.73 in Fig. EJ), and a running state, in which the particles surmount the 
barriers of the potential. Moreover, one can distinguish two classes of running states: 
either the particle overcomes the barriers without any back-turns (cf. the case with an 
amplitude a = 3.70 in Fig. |3J) or it undergoes frequent oscillations and back-scattering 
events (cf. the case with amplitude a = 4.55 in Fig. EJ). For a small driving amplitude, 
we find that the locked behavior is generic implying that the average motor velocity is 
almost zero, see Fig. If the amplitude is increased up to some critical value, here 
a = 3.25, the running solutions emerge. Around that critical point, there occurs a 
'battle of attractors' and the particle burns energy for both barrier crossings and intra- 
well oscillations. This behavior is reflected in an enormous enhancement of the effective 
diffusion jTHj. 

If the driving amplitude is further increased, a regime of optimal transport sets 
in. The rapid growth of the average velocity is accompanied by a decline of both the 
position and the velocity fluctuations. It means that the trajectories bundle closely 
together; note the case a = 3.70 in Fig. Because there are no intra-well oscillations, 
the energy that gets dissipated per unit distance, is minimal. 

At even larger drive amplitudes an upper threshold is approached (in the present 
case this threshold is located at around a = 4.7) where the velocity sharply decreases 
to a value close to zero. Moreover, the diffusion coefficient is small and the velocity 
fluctuations are large, cf. the case with amplitude a = 4.73 in Fig. 03 In this regime, 
the particle dangles around its actual position, as it occurs for a < 3, meaning that its 
motion is confined mostly to one well. We note, however, that the amplitude of the intra- 
well oscillations becomes much larger so that the corresponding velocity fluctuations are 
also large. 

We conclude that the diffusion coefficient is small for cases when the particle 
performs either locked motion or running motion without back-turns. 

All these considerations are accurately encoded and described by the two previously 
discussed measures, namely, the efficiency and the Peclet number Pe in © or in 
(114)1 . It is found that the optimal regime for the ideal modus operandi of the Brownian 
motor is achieved when both the efficiency and the Peclet number become maximal, see 
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Figure 5. (Color online) Brownian trajectories of the rocked particle moving in the 
asymmetric ratchet potential V(x) — Vq [sin(27nc) + 0.245 sin(47rizi) + 0.04 sin (67ra;)], 
where Vo — 0.461 normalizes the barrier height to unity. The forces stemming from 
such a potential range between —4.67 and 1.83. The two angular frequencies at the well- 
bottom and at the barrier-top are the same, reading 5.34. The remaining parameters 
are: 7 = 0.9, co = 4.9 and Do = 0.001. The values of the driving amplitude are 
a = 3.29,3.37,3.70,4.41,4.55,4.73. One can see that for a = 3.29 and 4.73 in part 
(a) (in blue online), the particles usually oscillate in a potential well, most of the time 
performing only a few steps. This results in an almost zero mean velocity, a very small 
effective diffusion but with rather large velocity fluctuations. For another set of driving 
amplitudes: a — 3.37 and 4.55 in part (b) (in green online) the mean velocity is large, 
o~ v becomes suppressed, but the effective diffusion exhibits an enlargement due to a 
'battle between attractors'. Part (c): The cases a = 3.70 and a = 4.41 (in red online) 
correspond to the optimal modus operandi of the inertial Brownian motor - the net 
drift is maximal and fluctuations get suppressed. 
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Figure 6. (Color online) Left top: the average, dimensionless velocity (v) of 
the inertial, rocked Brownian motor under nonadiabatic driving conditions. Left 
bottom: corresponding velocity fluctuations a v (dotted line) and corresponding 
diffusion coefficient D e ff (solid line). Right top: Brownian motor efficiency rj. Right 
bottom: depicted is the Peclet number Pe, being proportional the inverse of the Fano 
factor. All quantities are plotted versus the external driving amplitude a. Values of 
the remaining parameters are the same as in Fig. El e.g. the thermal noise strength 
here is Dq = O.OOf . The numerical errors are within the line width. 



in Fig. El Indeed, in this regime of optimal performance, the particle moves forward 
steadily, undergoing rare back-turns [TH], see the case a = 3.70 in Fig. |3 

6. Role of temperature 

We next address the dependence on the strength of thermal noise. In Fig. [?1 we present 
our numerical results for the noise-assisted, directed transport at a larger temperature, 
namely for D = 0.005. The potential barrier height is still rather high in comparison 
to the thermal energy In this regime, a so-called current reversal, i.e. a change of 
the transport direction occurs as a function of the driving amplitude. Otherwise, the 
behavior remains qualitatively the same as for lower noise D = 0.001. The diffusion 
coefficient exhibits three maxima and two minima in the corresponding interval of the 
drive amplitudes. However, the optimal regime corresponds to the neighborhood of the 
second minimum of the diffusion coefficient. In contrast, we notice that at lower noise 
D = 0.001, the optimal regime set in within the neighborhood of the first minimum of 

Deff- 

Under higher temperature operating conditions, optimal transport also occurs when 
both the efficiency 77 and the Peclet number Pe are maximal. 




7. Summary 

In this work, criteria for the optimal transport of an inertial rocked ratchet were 
established using two characteristic quantifiers: the energetic efficiency Q and the 
Peclet number (JBJ). Adapting the methods of point processes to rocked Brownian 
motors, we expressed the averaged motor velocity and the position-diffusion coefficient 
by corresponding averages of the point process Both these measures can be obtained 
from simulations of the driven Langevin dynamics (J2J). 

The Fano factor F used in the theory of point processes is related to the Peclet 
number in a simple manner via Pe= 2/F. In our case, it is more convenient to employ 
the Peclet number because in regimes where the average velocity is very small, the Peclet 
number assumes values close to zero, while the Fano factor would diverge. From our 
numerical analysis it follows that the optimal modus operandi for the inertial Brownian 
motor is obtained when the efficiency rj and the Peclet number assume simultaneously 
maximal values. 



Appendix 

In the present paper we have considered the effective diffusion coefficient, which is 
defined as 

where the brackets (...) denote an average over the initial conditions of position and 
velocity and over all realizations of the thermal noise. Another definition of the diffusion 
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coefficient is given by the formula 

D=]im <Mt)- Sx (o)?) i 

t^co 2t 

where 5x(t) = x(t) — (x(t)). By inspection one finds 

Deff = D (17) 

if 

lim -(5x(t)6x(0)) = 0. (18) 

t->00 t 

In our case, this term vanishes because of the presence of thermal noise and dissipation. 
More generally, \{5x(t)8x(0))\ may increase at most as t 1 ^ 2 if the diffusion coefficient 
D as defined in (JTHJ) is finite. Consequently, for such processes the equation (fT7|) also 
holds. 

We now show that the diffusion constant D is related to the auto-correlation 
function of the velocity via a Green-Kubo relation, in spite of the fact that the system 
is far from equilibrium. For a system with periodic driving, D takes the form 

D= / dsC(s), (19) 
Jo 

where 

C(s) = ^ ( T drC as (r,s) (20) 
1 Jo 

denotes the time average of the velocity correlation function C as (r,s) over one period 
T — 2-k/uj of the driving and where 

C as (t, s) = (5v(t)5v(t + s)) as = (5v(t + s)5v(t)) as (21) 

is the nonequilibrium asymptotic velocity-velocity correlation function. In the case of 
periodic driving, this function is periodic with respect to the first argument, i.e., 

C(t,s) = C(t + T,s). (22) 

To show the Green-Kubo relation, we start from the expression x(t) = v(t) from which 
it follows that 

Sx(t)-5x(0) = I ds6v(s). (23) 
Jo 

Therefore (fTBj) takes the form 

I rt rt 

2D = lim - / dsi / ds 2 (5v (si)Sv(s 2 )) 

t^ca f Jo JO 
1 rt rt 

= lim - / dsi / ds 2 C(s 2 , si - s 2 ), (24) 
t^oo t Jo Jo 

where 

C(t,s) = (6v(t)6v(t + 8)). (25) 
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Changing the integration variables (si, s 2 ) — > (s = Si — s 2 , r = s 2 ) and exploiting the 
symmetry of the correlation function, C(t, s) — C(t + s, — s), one obtains 

■■ [ S drC{r,s) 
Jo 

C(t, s) - lim - C ds f dr C(r, s). (26) 

t-^oo t Jo Jt-s 

We assume that the diffusion coefficient is finite. Therefore the second term in the 
second line of (J2T?|) tends to zero as t —>■ oo, so that 

POO 1 ft 

D= ds lim- / drC(r,s). (27) 
jo t Jo 

For t = i^T, one splits the second integral into sum over subsequent periods, 

I pt 1 K fkT 

lim - / drC(T,s)= lim — — / drCir^s) 
t^oof Jo K 1 k^oo KT ^ 7(fc-i)T v ; 



T Jo 



drC as (r,s) = C(s) (28) 



where 



C as (r, s ) = lim 1 ^ ^(^ + fcT ' s )- ( 29 ) 

The Eqs. ()25|) . ()27 p -()29 p represent the Green-Kubo relation for the diffusion constant 
of such periodically driven processes x(t); notably, these per se constitute far from 
equilibrium processes. 
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Abstract. The expression for the effective diffusion of an inertial, periodically driven 
Brownian particle in an asymmetric, periodic potential is compared with the step 
number diffusion which is extracted from the corresponding coarse grained hopping 
process specifying the number of covered spatial periods within each temporal period. 
The two expressions are typically different and involve the correlations between the 
number of hops. 



The expression used for the diffusion constant D c $ in eq. (12) in Ref. PQ, which will 
be denoted by in the sequel, coincides with the definition in eq. (5) in Ref. [I| only 
under certain conditions, see the Fig. |HJ In order to understand the relation between 
these expressions we split the random distance x(nT), which the particle has covered 
after n periods of duration T, into its integer multiple N(nT) of spatial periods of length 
L and a remainder e(nT) 

x(nT) = N(nT)L + e(nT) (El) 

Note that e(nT) is non-negative and bounded by L. From the definition (5) we then 
obtain 
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Each fluctuation SN(nT) = N(nT) — (N(nT)) contributes to the averages with a factor 
growing as n 1//2 such that only the first term on the right hand side of the first equation 
contributes in the limit n — > oo. Next, we represent the number N(nT) of spatial 
periods covered within n temporal periods T as the sum of the number iVj. of spatial 
periods which the particle passes through within the fcth temporal period of the driving 
force, i.e. 

N(nT) = J2N k (E3) 
k=i 



Optimal strategy for controlling transport in inertial Brownian motors 



2 



From eq. 



21) we then obtain 
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is the quantity that was used in eq. (12) in Ref. [JJ. Here we took into account that 
in the limit n — > oo the increments 5N k become stationary and the variances ((5Nk) 2 ) 
are independent of k. The difference between D e s and Dn results from the sum over 
the correlations between the increments 5N k . In the limit n — > oo the double sum is 
dominated by terms with large values of k and I for which the correlations (SNkSNi) 
only depend on the difference k — I. If the correlations decay faster than (k — /)" 2 the 
limit can be simplified to read 
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Figure 8. The step number diffusion Dm and the effective diffusion -D e // are depicted 
in panel (a) as functions of the driving amplitude a on a logarithmic scale for the driving 
frequency uj = 4.9, noise strength D — 0.001 and potential parameters c\ — 0.425 and 
c-2 = 0.04. In the locked regime for a < 3.1 both diffusion constants are comparably 
small. In the running regime for a > 3.1 both diffusion constants are large with D e g 
becoming larger than Dn roughly by a factor of ten. The normalized correlations 
(SNi +n SNi) I ((5Ni) 2 ) are depicted in panel (b) for two different values of a, which 
are marked by arrows in panel (a). They extend over many periods T leading to the 
observed discrepancy of the two diffusion constants. For a — 3.25 the decay of the 
correlations is much faster than for a = 3.7. Accordingly the difference D c ff — Dn is 
more pronounced at the larger a value. 
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In principle, this sum may take positive as well as negative values. In the Fig. |S] 
we display the dependence of D efi and D N on the amplitude of the driving force and 
illustrate illustrate the correlation functions (5Nk+ m SNk) in two selected cases. 
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